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ABSTRACT 

We explore several physical effects on the power spectrum of the Lya forest trans- 
mitted flux. The effects we investigate here are usually not part of hydrodynamic 
simulations and so need to be estimated separately. The most important effect is 
that of high column density absorbers with damping wings, which add power on large 
scales. We compute their effect using the observational constraints on their abundance 
as a function of column density. Ignoring their effect leads to an underestimation of 
the slope of the linear theory power spectrum. The second effect we investigate is that 
of fluctuations in the ionizing radiation field. For this purpose we use a very large 
high resolution N-body simulation, which allows us to simulate both the fluctuations 
in the ionizing radiation and the small scale Lya forest within the same simulation. 
We find an enhancement of power on large scales for quasars and a suppression for 
galaxies. The strength of the effect rapidly increases with increasing redshift, allow- 
ing it to be uniquely identified in cases where it is significant. We develop templates 
which can be used to search for this effect as a function of quasar lifetime, quasar 
luminosity function, and attenuation length. Finally, we explore the effects of galactic 
winds using hydrodynamic simulations. We find the wind effects on the Lya forest 
power spectrum to be be degenerate with parameters related to the temperature of 
the gas that are already marginalized over in cosmological fits. While more work is 
needed to conclusively exclude all possible systematic errors, our results suggest that, 
in the context of data analysis procedures where parameters of the Lya forest model 
are properly marginalized over, the flux power spectrum is a reliable tracer of cosmo- 
logical information. 



1 INTRODUCTION 



One of the primary goals of cosmology is the determination 
of the origin of structure in the universe. Among the avail- 
able probes today one of the most promising is the amplitude 
and shape of the power spectrum of the primordial fluctu- 
ations. It can be constrained both by the CMB and other 
probes of large scale structure. Combining different data sets 
that cover a large range in scale is particularly powerful. Of 
the current cosmological probes, the Ly-a forest - the ab- 
sorption observed in quasar spectra by neutral hydrogen in 
the intergalactic medium (hereafter IGM) - has the poten- 
tial to give the most precise information on small scales. 
It probes fluctuations around megaparsec scales at redshifts 
between 2-4, so nonlinear evolution, while not negligible, has 
not erased all of the primordial information. The statistic of 
choice when analyzing the Lya forest is t he power spectrum 
of the transmitted flux fraction, Pf(&, z) iCroft et al.ll99slK 
While several groups have extracted this quanti t y from the 
data JCroft et al.ll2002bl: iMcDonald et alJl2000t Kim et alJ 



2003), a recent analysis of 3300 quasar spectra from the 
Sloan Digital Sky Survey (SDSS) has i mproved the sta- 
tistica l errors by an order of magnitude jMcDonald et all 
2004a). This significant increase in statistical power must be 
accompanied by a corresponding improvement in the theo- 
retical calculations if the new data are to be exploited to 
the fullest extent. The goal of this paper is to make a step 
in this direction. We focus on three topics: high density sys- 
tems (especially the effect of damping wings), fluctuations in 
the ionizing background, and effects of galactic superwinds. 

In the standard picture of the Ly-a forest the gas in the 
IGM is in ionization equilibrium. The rate of ionization by 
the UV background balances the rate of recombination of 
protons and electrons. The recombination rate depends on 
the temperature of the gas, which is a function of the gas 
density. The temperature-density relation can be parameter- 
ized by an amplitude, To, and a slope 7 — 1 = dlnT/dlnp. 
The uncertainties in the mean intensity of the UV back- 
ground, the mean baryon density, and other parameters that 
set the normalization of the relation between optical depth 
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and density can be combined into one parameter: the mean 
transmitted flux, F(z). An additional nuisan ce parameter 
is the filtering length fcp llGnedin fc Hui 1998), which is re- 
lated to the thermal history of the IGM. The parameters of 
the gas model, Jo, 7 — 1, F, and ftp, are marginalize d over 
when c omputing constraints on cosmology llMcDonald et all 
l2004bh . Thus any additional physical effects must be differ- 
ent from the combined effect of these 4 parameters if they 
are to be relevant for the cosmological analysis. Our analysis 
differs from previous analyses in that we include this factor 
when assessing the importance of a given physical effect. 

The first effect we investigate is that of high density 
absorbers on the flux power spectrum. In the bulk of the 
intergalactic medium (IGM) at high redshift, neutral gas is 
in ionization equilibrium, with the proton-electron recom- 
bination rate balancing ionization of neutral hydrogen by 
a fairly uniform ionizing background. Systems with neutral 
column density A(Hi ) > few x 10 17 atoms cm -2 become 
self-shielded, i.e., the exterior of the system absorbs all in- 
coming ionizing radiation so the interior no longer sees the 
background, becoming mostly neutral. These absorbers are 
called Lyman-limit s ystems (LLS). LLSs should be relatively 
small (Schave 2001), and for that or other reasons, such 
as too small simulation boxes, are no t generally well repro 



duced in hydrodynamic s i mulations l|Miralda-Escude_et_al 
| 1996|: iGardner et aHfeOOlt ICen et al.lf2003t iNaeamine et al 
2004). As the column density of the systems increases to- 
ward the traditional definition of a dam ped-Lyq absorber 
(DLA), AT (Hi ) > 2 x 10 20 atoms cm -2 iwblfe et alJll986t 
ISmith et alJll986ft , damping wings come to dominate the 
equivalent width of the systems, enhancing their impact on 
spectra. In addition to the problem of simulations not nec- 
essarily reproducing the number of systems, damping wings 
are not usually included in the simulated spectra used to 
predict Pp(k,z) at all, because the widest of them can ex- 
tend to the full width of a typical simulation box. 

The effect of damping wings differs from the effects 
of UV background fluctuations and galactic winds, dis- 
cussed next, in that they are certainly present in the data, 
and the systems are even more or less directly observ- 
able, although the detailed column density distribution be- 
low 2 x 10 20 atoms cm -2 is more difficult to resolve ob- 
servationally and has not been the subject of much in- 
vestigatio n. The poss i ble im portance of DLAs was recog- 
nized by ICroft et al] il999ft , who investigated their effect 
by measuring dPp/dk with and without the identified sys- 
tems in their spectra, finding a change t hat was generally 
not larger than their statistical error bars. IViel et all ]2004f) 
emphasized the non-negligible contribution of systems with 
AT(H i ) > 10 15 atoms cm -2 to the flux power. In this paper 
we investigate the effect of high density systems in general in 
simulations, and particularly the importance of the damping 
wings of absorbers with A(Hi ) < 2 x 10 20 atoms cm" 2 . 

The second effect we analyze in this paper are the 
fluctuations in the ionizing background. The intensity of 
the ionizing background determines the density of neu- 
tral hydrogen and if the intensity is spatially varying this 
will lead to spatial variations in neutral hydrogen den- 
sity. These intensity fluctuations in the ionizing background 



the past JZudll99l 


Fardal & Shullll993l ICroft et alJll999: 


Gncdin & Hamilton 


2002). The two most recent works in- 



vestigated it in some detail focusin g ^ on Lya forest statistics 
such as the flux power spectrum IIMeiksin fc White! 120041 : 
ICroftll2003l) . iMeiksin fc White! J20041) use relatively small 
dark matter PM simulations to investigate the effect, focus- 
ing on high redshifts (z > 4) where attenuation lengths are 
short and the small simulation box sizes used are adequate. 
They assume quasar hosts are uncorrelated with the Lya 
forest and as a consequence they find the effe ct to enhance 
the fluctuations on large scales. In contrast, ICroftl j2003fl 
used large P 3 M dark matter simulations in combination 
with small scale hydro-dynamic simulations to investigate 
the effect at z = 3. In this model quasar hosts are placed 
in high density regions and are correlated with the Lya for- 
est. Since there is a higher intensity ionizing background in 
the regions where there is also more neutral hydrogen, a 
cancellation occurs and fl uctuati o ns in the Lya forest are 
suppressed on large scales. ICrofj (2003) introduced several 
additional improvements in the modeling of the effect, such 
as inclusion of light cone effects due to quasar finite lifetime, 
photon shadowing by high density neutral hydrogen along 
the line of sight, which was investigated using a photon ray 
tracing code, and an investigation of the effects of beam- 
ing. While the photon shadowing and beaming appear to 
be of little relevance for the Lya forest flux power spectrum 
statistics, the light cone effects increase the effect by up to 
a factor of two. In this paper we will therefore ignore shad- 
owing and beaming, using the uniform isotropic attenuation 
approximation, but we will include the lightcone effects. 

The main numerical issue when investigating the effect 
of ionizing background fluctuations on the Lya forest statis- 
tics is the required dynamic range of simulations. The ioniz- 
ing background attenuation length in comoving units ranges 
between 50 h~ x Mpc and 500 h~ x Mpc in comoving units, so 
to properly estimate the fluctuations box sizes of order sev- 
eral hundred megaparsecs are required. On the other hand, 
to properly simulate the Lya forest one needs resolution be- 
low 100fo -1 kpc to resolve the Jeans scale. Combining the 
two requirements leads to a dynamic ra nge not ava ilable 
presently in hydrodynamic simulations. In lCroftJ i2003t) this 
problem was approached using a hybrid method where large 
scales are simulated with a dark matter only simulation and 
a hydrodynamic simulation is used to simulate small scales, 
by randomly choosing a patch that matches the P 3 M simu- 
lation in the value of the density field smoothed at a larger 
scale. However, this approach does not preserve all of the 
correlations present in a single full resolution simulation and 
it is not clear if the results can be used for quantitative es- 
timates. 

In this paper we use pure dark matter simulations 
to simulate both ionizing background fluctuations and the 
Lya forest. We perform the analysis as a function of red- 
shift over the range 2 < z < 4.6. Dark matter simula- 
tions ca n reproduce qualitat i vely full hydrodynamic simu- 
lations llGnedin fc Huilll998l; IMeiksin fc White! [200 ll) . Our 
approach is similar to IMeiksin fc White! i2004ft . where it 
was shown, using a series of resolution studies, that to as- 
sess the relative effects between simulations with and with- 
out fluctuations it suffices to have a 128 3 simulation in a 
30 ft," 1 Mpc box, even if the small scale fluctuations in the 
Lya forest are under-resolved in an absolute sense. In this 
paper wc use a 1024 3 TPM simulation in a 320 h' 1 Mpc box 
bode fc Ostrikedl2003ft . whose mass resolution and force 
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resolution in low density regions is similar to a 128 3 PM 
simulation in a 30 h -1 Mpc box, while its force resolution in 
high density regions is significantly better. The size of the 
simulation is several attenuation lengths at z — 3 and we use 
periodic boundary conditions to extend the calculations of 
photon propagation from quasars to even larger distances. 
Thus we expect this simulation to be sufficiently accurate 
for our purposes both on large and small scales. An addi- 
tional advantage of its high mass and force resolution is that 
it resolves halos with masses above W 11 h' 1 M Q . This allows 
us to generate a halo catalog using a friends-of-friends halo 
finder, which we used as a catalog of quasar hosts. 

The third physical effect on the Lya forest that we ad- 
dress in this paper is that of galactic superwinds (GSW). 
There is clear evidence that powerful winds emerging out o f 
galaxies are present bo th at low re dshifts jHeckmanlbOfjil 
and at high redshifts ijPettini et all 1200 ilk but the extent 
to which these affect the intergalactic medium (IGM) at 
large separations from the host galaxy is poorly known. Cur- 
rent studies of Lya forest absorption in background quasars 
at small angular separations from a Lyman-break galaxy 
(LBG) suggest that winds (or something l ike them) may be 
affect ing IGM to an appreciable distance l|Adelberger et all 

2003) , but the statistics are poor and there could be se- 
lection effects involved. However, there is also ample ev- 
idence of metals in the i ntermediate to low density IGM 
JSongaila fc Cowid Il99fit ISchave et all 1200.4 Uracil et all 

2004) , which suggests that GSW must be operating in trans- 
porting metals from galaxies to the lower density IGM. 

There have been several recent sim ulations inve stigat- 
ing the effect of GSW on the Lva forest iCroft et alJl2002eJ: 
[ McDonald et al]l2002l iBruscoli et al]l2003f Kollmeier et alJ 

1200.4 r 



2002, _ 

1200.4 iMaselli et alJl2004HDesiacaues et al.ll200.4 . Many of 
these focused on the correlations with L BGs, attempting 
to rep roduce the observational results of lAdelberger et alJ 

(2003) . While they are generally successful in reproducing 
the observations at larger distances, they fail in the inner 
~ /i _1 Mpc, where the observations suggest that there is 
frequently little or no absorption, while simulations would 
suggest that the absorption is significantly larger because 
of the inc reased dens i ty of n e utral hydrogen close to the 
galaxy. In lCroft et all J2002aft : IPesiacaues et all fcOOSl) the 
effects of winds on the Lya forest flux p o wer sp ectrum have 
also been investigated. For ICroft et alJ <l2002al) the effects 
are small, but not negligi ble at the level require d by new 
Lya forest measuremen ts (IMcDonald et al1l2004al) . while in 
iDesiacaues et all (12003ft the effects appear to be completely 
negligible. To explain the observations the winds were as- 
sumed to blow out a bubble around LBGs, where the bubble 
radius was put in by hand. Some of these comparisons were 
based simply on comparing wind and no wind simulations to 
each other, while winds produce other effects such as heating 
of the IGM. In this case it is not clear if the effects produced 
by the winds are unique or degenerate with changes in the 
mean absorption level, temperature-density relation, or fil- 
tering length, all of which we marginalize over in the stan- 
dard analysis anyways. In this paper we will address these 
issues by looking for residual effects after these effects are 
accou nted for. We will use wind simulations from lCen et alJ 

(2004) , which have been calibrated to the observed wind 
velocities, but more extreme cases are also explored. 



2 HIGH DENSITY ABSORBERS 

In this section we investigate the effect of high dens i ty ab - 
sorb ers on Pf(fc , z). T he investigations bv lCroft et alJ (Il999l) 
and IViel et alJ <2004) were not performed at the level of 
de tail we need to c ompar e to t he P_F(fc,z) meas u rement 
of iMcDonald et alJ ll2004al) . The iMcDonald et alJ (l2004al) 
measurement an d the ac companying simulation analysis by 
McDonald ct al spans the redshift range 2.2 ~< 

z ~< 4.2 in bins spaced by Az = 0.2, and probes scales 
0.0013 s/km < k < 0.02 s/km in bins spaced by A log 10 (fc) = 
0.1. The precision of the 132 resulting Pf(&, z) points is usu- 
ally better than 10% and often as good as 3%. Our goal in 
this section, once we find an effect significant enough to re- 
quire us to account for it when fitting to the data, is to 
create templates for the dependence of the effect on k and 
z, with a level of accuracy sufficient to match the present 
data. 

Within the context of the model for the IGM that we are 
assuming, as represented in the hydrodynamic simulations, 
systems with JV(Hi ) < 1.6 x 10 17 atoms cm" 2 should not 
present any problem for us, so we will focus on higher col- 
umn density systems. The iV(Hi ) < 1 .6 x 10 17 atom s cm -2 
systems are close enough to resolved dSchavell200i in sim- 
ulations that if they are important to the power spectrum, 
in a way that is not yet fully resolved, the power spectrum 
would change as the resoluti on is increased in s tandard tests 
such as those performed in IMcDonald et alJ ll2004bl) (i.e., 
there is no reason not to expect a smooth change in the 
power spectrum with increasing resolution). Systems with 
iV(Hi ) > 1.6 x 10 17 atoms cm -2 seem more likely to be a 
problem because the onset of self-shielding causes a rather 
sudden increase in optical depth with density, the impor- 
tance of which might not be detected by resolution tests 
(this is not obvious, but is conceivable). 

The section is arranged as follows: In H2.ll we explain 
the observationally determined column density distribution 
that we use in H2.2l and l2.3l Then in t|2.2l we look at high den- 
sity absorbers in the hydro dynamic simulations used to pre- 
dict the power spectrum in IMcDonald et all (l2004bl) . with- 
out considering damping wings. Finally, in H2.3I we estimate 
the effect of the damping wings of the highest density sys- 
tems on Pp(k, z). 



2.1 Estimate of the log(TV) > 17.2 Column Density 
Distribution from Observations 

Subsections 12.21 and 12.31 rely on an estimate of the col- 
umn density distribution, f(N), for all iV(Hi ) > 1.6 x 
10 17 atoms cm -2 . Traditionally, the full distribution has 
only been measured for iV(Hi ) > 2 x 10 20 atoms cm -2 , 
with the LLSs grouped together into a total number density 
n(z), because the column densities of these systems are hard 
to measure. In the course of our investigation of damping 
wings, we determined that the systems with log(TV) < 20.3 
were significant, so we decided to re-estimate the full column 
density distribution from the data. 

The number of systems as a function of column den- 
sity, f(N), is expected to have a non-trivial form in the 
LLS-DLA regime, because of self-shielding. If one imagines 
increasing the density of a reference system, the implied 
neutral density increases gradually in ionization equilib- 
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rium with the background radiation until the system reaches 
log(iV) ~ 17.2, when self-shielding reduces the effective ion- 
izing background and the neutral density suddenly increases 
much more ra pidly until the system becomes almost com- 
pletely neutral. Zheng & Miralda-Escude ( 2002) give predic- 
tions for the resulting column density distribution, based on 
a spherical isothermal halo model. The model has two free 
parameters: the overall normalization of f(N) (proportional 
to the number density of halos), and a parameter that we 
will nominally call the mass of the halos, although this in- 
terpretation should not be taken literally as the value of the 
parameter depends o n other things like the strength of the 
ionizing background Iph eng fc Miralda-Escudel I2002L The 
IZheng fc Miralda-Escudel (I2002T) model is a toy model, but 
it seems reasonable to assume that the class of f(N) shapes 
that it describes contains a shape close enough to the truth 
for our purposes. 

The two parameters, amplitude and halo mass, are not 
predicted by the model so we will determine them by a fit 
to the data. To describe the column density distribution 
including dependence on redshift, we perform a maximum 
likelihood fit using the formula 

f(N,z) = A xl A f XK (N,M* x™) , (1) 

where x z = (1 + z)/(\ + z+), and fzM is the fu nctional 
form described bv lZheng fc Miralda-Escudel J2002t) . i.e., we 
allow the overall normalization to evolve as a power law 
in 1 + z, and allow the mass of the typical halo to evolve 
similarly. In practice, we include the M dependence by 
simply interpolating between t he four example curves in 
IZheng fc Miralda-Escudtl ll2002h 's Figure 2 (using linear in- 
terpolation in ln(/) and ln(iV)), which were provided in nu- 
merical form by Z. Zheng. 

We take the data in the DLA regime from 
Prochaska & Hcrber^Tbrn 12004]). For reference, 
Prochaska fc Herbert-Forti 120041) estimate the num- 
ber density of DLAs per unit redshift, finding 
riDLA(z) = 0.189 ± 0.026 for a redshift bin with 
2.0 < z < 2.5, 0.215 ± 0.032 for z = 2.5 - 3.0, 0.271 ± 0.049 
for z = 3.0 - 3.5, 0.366 ± 0.073 for z = 3.5 - 4.0, and 
0.4 01 ±0.121 for z = 4 — 5. We d o not use these numbers, or 
the iProchaska fc Herbert-Forti i2004T) column density dis- 
tribution directly, because they would not give an optimal 
constraint on our model parameters. Instead, our maximum 
likelihood analysis uses the column densities of individual 
absorbers and the segments of quasar spectra searched, 
considering only z > 2.0 (these data were provided in 
machine reada ble form by J. Prochaska). 

Recently, IPeroux et alj l2003al) made a first mea- 
surement of the abundance of systems with 1.0 x 
10 19 atoms cm~ 2 < N(Ri ) < 2.0 x 10 20 atoms cm" 2 , find- 
ing n 19 -2o.'i(z < 3.5) = 0.49±? and 7219-20.3(2: > 3.5) = 
0.97±? (there were a total of 7 and 3 systems, respectively, 
from which one can estimate the unsupplied error bars) . We 
use their data quasar- by- quasar, reading from their Table 1, 
in the same way that we use the DLA data. 

Finally, we add a constraint on the integra l of f(N) over 
the LL S column density range, tills (z) from IPeroux et 
(2003b). They measured nLLs(.z) = no(l + z)"' LLS per unit 
redshift, with n LLS = 0.07±g ol and 7 = 2.45±g;e|, for LLS 
with 2.4 ~< z ~< 4.5; however, note that the error on no 
is useless because the correlation with 7 is not given - the 



error on n(z) at the center of weight of the data will be 
smaller than that implied by the errors on no ■ IPeroux et alJ 
l2003tft 's sample included 67 LLSs, implying that the error 
on the number density should be about 12%. 

The results of the maximum likelihood fit are A = 
25.1 ± 3.6, 7 A = 1.80 ± 0.80, log(M/M Q ) = 9.36 ± 0.18, 
and 7aj = 1.4 ± 2.4, for z* = 3.0. We give the best fit 
numbers only so that the column density distribution we 
use will be reproducible. Our application only requires that 
something resembling the true column density distribution 
falls within the class of curves explored by varying M, which 
seems likely - the interpretation of M as a halo mass is much 
more complicated, as discussed bv lZheng fc Miralda-Escudel 
( 2002) (among other problems, these curves were computed 
for an Einstein-de Sitter Universe). Note also that the er- 
rors are generally correlated. Performing an analysis of this 
type aimed at better understanding LLSs and DLAs through 
their column density distribution would be an interesting 
topic for future work. 

2.2 Self-shielding and High Density Absorbers in 
Simulations 

We start by examining the case where we ignore damping 
wings, because the wings require a different treatment using 
mock spectra that can be much longer than the width of 
the simulations. We will see that this is ultimately an aca- 
demic exercise, because the damping wings turn out to be 
the most important aspect of high density systems, but it is 
informative nonetheless. 

Th e hydrodynamic simulations used in lMcDonald et all 
(2004b), which we will focus on in this subsection, are de- 
scribed in more detail in lCen et al.H2003l) . We simulate a flat 
ACDM cosmological model with Q m = 0.3, Qb = 0.04, and 
h = 0.7, using an L = 10 hT x Mpc box, with the gas proper- 
ties tracked in N — 256 3 fixed grid cells and the dark matter 
traced by 256 3 particles. While lCen et alJ 12003T) reproduces 
the DLA abundance after allowing for dust obscuration, very 
high density systems in general are not necessarily well re- 
produced by ours or any other hydrodynamic simulations 
ICen et all2003t iMiralda-Escude et alJll99rilGardner et alJ 
l200ltlNagamine et al.l2004D . Even if they were, we could not 
compute the effect of strong damping wings on our power 
spectrum accurately because of the limited width of our sim- 
ulation boxes. 

We start by estimating the effect of self-shielding in our 
hydrodynamic simulations. The simulations were performed 
with a rough self-shielding approximation: the background 
radiation for each cell was attenuated by the column density 
of that cell. Assuming ionization equilibrium, we can recom- 
pute the neutral densities with and without this shielding 
(in the shielded case three iterations of the density calcu- 
lation, starting with no shielding, are required for conver- 
gence) . The effect of self-shielding on the amount of neutral 
gas at a = 0.24 is dramatic - an increase by a factor 70; 
however, in the absence of damping wings the effect on F 
is tiny - 0.67561 is reduced to 0.67559. The differences are 
the same size (~0.00002) at a = 0.2 and 0.32, and the dif- 
ference in Pp(k,z) is never more than 0.1% over the range 
0.0013 s/km < k < 0.05 s/km. This remarkable result is 
not hard to understand: in the absence of damping wings, 
the equivalent width of an absorption line susceptible to self 
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shielding, i.e., with N(H.i ) > 1.6 x 10 17 atoms cm -2 , in- 
creases by only ~10% per order of magnitude increase in 
AT (Hi ), and these lines are very rare to begin with. This 
is an example of a more general fact about the Lya forest 
- it simply is not very sensitive to the details of the high 
density regions of the Universe, except, as we will see, when 
damping wings become important. 

Our simulations in fact produce too few LLS, even with 
self-shielding. To be sure that the missing systems do not 
significantly affect our results, we tried simply arbitrarily 
increasing the column density of the most dense systems 
in the simulation so that we match the observed number 
of LLS. The affect is slightly larger than the affect of self- 
shielding, but still tiny - F decreases by 0.0004, and Pf(A;, z) 
changes by less than 0.8%. 

In summary: Without damping wings, LLS and denser 
systems are completely irrelevant to the use of Pf(&,z) for 
cosmology. Self-shielding, despite changing the amount of 
neutral gas by orders of magnitude, does nothing to the bulk 
of the IGM, and there is no significant different between the 
mean absorption or Pp(fe, z) predicted by our unmodified 
simulations and simulations where we have artificially in- 
serted the observed number of LLS. 

2.3 Damping Wings 

We now turn to the effect of damping wings. Tradition- 
ally, LLS have been defined to be systems with 1.6 x 
10 17 atoms cm" 2 < AT (Hi ) < 2.0 x 10 20 atoms cm" 2 , 
where AT (Hi ) is the neutral hydrogen column density, while 
DLAs have AT (Hi ) > 2.0 x 10 20 atoms cm" 2 . This cutoff is 
somewhat arbitrary since systems with AT(H i ) < 2 x 10 20 
still show substantial damping wings. For example, the 
equivalent widths of cold systems with log[A r (Hi )] =(17, 
18, 19, 20, 21) are (54, 171, 539, 1705, 5392) kms" 1 , while 
the widths at the same densities would be (164, 185, 204, 
222, 237) kms" 1 if they were purely thermally broadened at 
temperature ~ 50000 K. We will explore the effect of ab- 
sorption by these systems on Pf(&, z) and F, using the ob- 
servation constraints on their column density from the next 
subsection, combined with several models for their spatial 
distribution. 

The power spectrum of a completely random distribu- 
tion of systems is shown in Figure along with the power 
from DLAs only (the column density distribution used is de- 
scribed in >I2.1I| . We see that the effect is quite significant, 
and that more than half of it (on relevant scales) comes from 
the log(Af) < 20.3 systems. 

The region of a simulated spectrum that is being re- 
placed by the high density system may be relevant, as we 
saw in the case with no damping wings, We include this ef- 
fect using the lognorma l model mock spectra described in 
iMcDonald et al.N2004aTl . These are not perfectly realistic, 
but they have the advantage over full numerical simulations 
that they can be of arbitrary length so the full extent of 
damping wings can be included. Figure [5] shows the power 
added to the mock spectra when high density systems with 
the observed column density distribution are added to the 
mock spectra, either at random locations, or using a map- 
ping of the highest column density systems to the highest 
density maxima in the lognormal field. We see that the loca- 
tion of the systems does matter. The change in high-fc power 
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Figure 1. Pp(k, z) from a random distribution of LLS and DLA 
(upper lines) or just DLA (lower lines), each with the observed 
column density di stribution, relative to t he observed Lya for- 
est P F (k,z) from iMcDonald et al.ll2004'ah . Red (dashed), black 
(solid), and green (dotted) lines show z = 2.2, 3.2, and 4.2. The 
error bars indicate the fractional error on the observed Pp(k, z) at 
z = 3.2 (the errors at z = 2.2 are very similar, while the errors at 
z = 4.2 are ~ 2 times bigger). The vertical dott ed line shows the 
largest scale used for constraining cosmology by ( McDonald et al. 
l2004bl) . 



is removed when the systems are inserted at already high- 
density locations, presumably because this power is related 
to relatively low column density systems that are not signif- 
icantly broadened by damping wings and thus produce ab- 
sorption that is insensitive to column density, as seen above. 
In fact, the difference between the added power in the two 
cases i s not significant to cosmological fits (IMcDonald et all 
l2004bt) . apparently because the fractional effect is is rela- 
tively flat (i.e., ^-independent) and thus mostly degener- 
ate with F inste ad of the primordial pow er spectrum. The 
agreement in the lMcDonald et al.l J2004bl) fits between these 
two cases makes us confident that detailed, realistic simula- 
tions of high density systems will not produce a significantly 
different result. The primary effect comes simply from the 
self-correlation across the width of individual systems with 
damping wings, which is directly computable from the ob- 
served column density distribution. 

How should we use these results to correct for their 
influence on the determination of the linear theory power 
spectrum, Pt(fc)? The effect is essentially additive, so the 
form P t hcory(fc, z) = P sim (k,z) +Pdam P (fc, z) is most natu- 
ral, where Ptheory(fc, z) is the final prediction to be compared 
with observations, P s i m (fc, z) is the prediction by simulations 
without damping wings, and Pdam P (fc, z) is the power shown 
in Figure Pdamp carries some uncertainty of course. As 
discussed in q2.ll the statistical errors on the number of 
LLSs and DLAs are probably about 15%, and the statistical 
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Figure 2. Change in Pp(k, z) when LLS and DLA with the ob- 
served column density distribution are inserted into mock spectra 
of the Lye? forest, relative to the observed Lye? forest Pp(k,z). 
The upper curves show the case where the high density systems 
are inserted randomly (different from Fig. in that there was no 
background forest in that figure), while for the lower curves the 
LLS and DLA were inserted at the highest density maxima in the 
mock Lye? forest spectra. The line types and error bars are as in 
Fig. Q Note that a consistent systematic error that is la for any 
single point is very significant, because we have many z and k 
bins to average over. 



error on Pdamp should be proportional to this. However, be- 
cause most observations do not distinguish column densities 
of LLSs, the important range of column densities just below 
the DLA cutoff is not directly constrained to the 15% level. 
For this reason, along with general caution to allow for any 
errors in our calculation, we suggest applying a 30% error to 
the amplitude of Pdam P (fc, z) in fits, i.e., adding APdnm P {k, z) 
to the simulated power, and [(A - l)/0.3] 2 to X 2 - An even 
more cautious analysis might allow for redshift dependence 
of A; however, there is not any reaso n to think a break from 
the lZheng fe Miralda-Esc udc (2002) column density distri- 
bution template used in 42. II should be significantly redshift 
dependent, and one might even argue that the constancy 
of the ratio of of Pdam P (fc»z ) to Pobservedffc, z) (Fig. |H is 
unlikely to be a coincidence. IMcDonald et al.l d2004bl) find 
a value of A = 1.01 ± 0.24, i.e., the power spectrum fit 
not only prefers our estimated amplitude, but is beginning 
to constrain the amplitude to better than our conservative 
30% external constraint. One might think that a more care- 
ful study could reduce this error and thus reduce the errors 
on the final Pl(&) measurement; however. IMcDonald et ail 
(2004b) shows that improving the error to 10% does not no- 
ticeably improve the Pl (k) results (this may change as other 
parts of the Ph{k) analysis are improved). 



3 FLUCTUATIONS IN THE IONIZING 
BACKGROUND 

The two most recent studies of fluctuations in the ioniz- 
ing bar4cgjTOnid_j ; ndtheir_effect on the flux po wer spectrum 
are lMeiksin fc White! ll2004l) andlCrofij d2003l) . The two dif- 
fer in details: while iMeiksin fc White! (12004ft use an N -body 
(PM) simulation and are limited to small box sizes, ICrofj 

(2003) uses a larger P Z M simulation and patches onto it 
sections of a high resolution hydrodynamic simulation. A 
large simulation box is needed to explore fluctuations in the 
ionizing background from quasars, whose typical separation 
is tens of megaparsecs. High resolution is needed to resolve 
the Lya forest. Here we use a very large TPM simulation 
with 320/i~ 1 Mpc box size and 1024 3 resolution to obtain 
both a l arge box and a high reso lution within the same sim- 
ulation dBode fc Ost riker 2003). This allows us to explore 
the background fluctuations on large scales and gives cor- 
rect correlations between large and small scales. While the 
simulation does not converge in an absolute sense on small 
scales, it probably converges for the purpose of relative com- 
parisons between simulation r esults with and withou t ion- 
izing background fluctuations dMeiksin fc Whitel2004ll . The 
cosmological parameters of the simulation are f2 m = 0.3, 
H = 70km/s/Mpc and cr 8 = 0.93. 

To generate a realization of the quas ar background we 
first s elect the sites of quasar hosts. In IMeiksin fc White! 

(2004) positions of these were chosen randomly, but this ig- 
nores the possibility that th e Lyq forest and background 
fluctuations are correlated. In Croft ( 2003) quasar host sites 
were chosen based on overdensity criteria, because the res- 
olution of the simulation was insufficient to resolve galactic 
size halos. We have actual halo catalogs from a halo finder, 
which are reliable down to masses around 1O 11 /i" 1 M (the 
particle mass is 3 x 10 9 h" 1 Mq). We use these halos as po- 
tential quasar hosts. We assume all quasars have the same 
lifetime toso, whic h we will vary between 10 7 — 10 s years 
dSteidel et al JI2002I) . 

We model the quasar luminosity function as a double 
power law, 

$(L) = (L/L.)0 + (L/L.)* ' (2) 

with the faint end slope Qi = 1.58. For the bright end we 
explore two possibilities, 3 = 3 .41 as suggested by quasars 
with z < 2.3 l)Bovle et alJl200(t. and 3 = 2.58 fo und in high 
redshift quasars with z > 4 dFan etai1l200ll) . The latter 
choice is flatter, has more bright quasar sources and leads to 
stronger fluctuations. The B band magnitude corresponding 
to L, in equation|5|is assumed to be given by M* = —22.65 — 
2.5(1.36,z - 0.27z 2 ), wh ile = 0.36 x 10" 6 Mpc~ 3 mag _1 
dMeiksin fc Whitdl2003l) . 

To generate a realization of quasars in the simulation at 
a given redshift we take the quasar luminosity function and 
multiply it by the ratio of the age of the universe at that 
redshift, ^Hubble, to the quasar lifetime, tQso- This gives us 
the density of quasar hosts, each of which will be a host to an 
active quasar at any time with a probability of tQso/^Hubbie- 
We then assign halo hosts to quasar hosts assuming a mono- 
tonic relation between the quasar luminosity (when active) 
and halo mass. We generate a turn-on time i on for each 
quasar choosing randomly between and tHubbic- In general 
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we choose quasars down to Mb = — 22, which is fainter than 
the turnover magnitude M*. 

To compute the background radiation at any given posi- 
tion we add up contributions from all quasars. We do this by 
first computing the time of observation i b s for a given posi- 
tion along the line of sight. We compute the distance r from 
that position to each quasar and convert it to light prop- 
agation time iprop- To count the contribution from a given 
quasar one must have ton + t pTop < t ohB < t on + t prop + t qso , 
so that the radiation from the quasar is passing through the 
given position at the time of observation. T his procedur e 
accounts for the light-cone effects discussed in lCrofd J2003). 

Each quasar adds a contribution to the radiation field T 
proportional to L r~ 2 exp(— r/r a tt)- The attenuation length 
r att is determined by the overall level of absorption of the ra- 
diation by neutral hydrogen. The actual value is somewhat 
uncertain and there are various estimates in the literature 
jFardal fc Shu3ll993t lHaardt fc Madal Il99fj: ICrofd l2003t 
iMeiksin fc Whitei2004) ~ For the low density Lya forest sim- 
ulations can be used. In addition there is also a significant 
contribution from high column density damped and Lyman- 
limit systems, which are typically not properly reproduced in 
simulations and their contribution has to be estimated sep- 
arately. Here we will use r att = 1.7 x 10 4 (1 + z) -3 - 2 ft _1 Mpc 
iMeiksin fc Whitd l2004'). but we will also explore the more 
pessimistic scenario where the attenuation lengths are 50% 
shorter (at the lower limit of the various estimates in the 
literature). In general, shorter attenuation lengths lead to 
larger fluctuations, since nearby sources account for a larger 
fraction of the radiation in this case. 

The attenuation length increases rapidly with decreas- 
ing redshift and at low redshifts it can exceed the simulation 
box size. This means that considering only the contributions 
from within the box will underestimate the background and 
so overestimate the fluctuations (since we always normal- 
ize the background relative to the required level needed to 
achieve a given mean level of absorption). We solve this by 
stacking together several boxes, summing over the quasars 
in all of them. These additional boxes are identical real- 
izations of the central box and we use periodic boundary 
conditions to patch them together. In total we perform the 
summation over 27 boxes. Note that while the additional 
boxes are identical to the original box, the actual quasars 
contributing to a given position are not repeated, since the 
quasar lifetime is much shorter than the light travel time 
through the box and so no single quasar can contribute to a 
given point more than once. Examples of the radiation field 
along a line of sight are shown in figure [3] This example is 
for z = 3, 4qso = 10 8 yr and r a tt = 161/i _1 Mpc. We see that 
most lines of sight are very smooth with radiation density 
close to the mean, but every now and then the line of sight 
passes very close to an active quasar and the radiation den- 
sity increases. In the most prominent cases one can see the 
light cone effects and times of quasar turn-on and turn-off. 

To generate Lya forest spectra from a pure N-body sim- 
ulation we choose a line of sight and use cloud-in-cell inter- 
polation of the nearby particles onto a one dimensional grid 
consisting of 1024 cells of the same size as the original 3- 
d simulation cells. This gives us the density and velocity at 
each grid cell. We compute the neutral density /optical depth 
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Figure 3. Examples of radiation background fluctuations along 
a line of sight assuming all of the background comes from quasars. 
This example is for z = 3.29 with r a tt = 161/i -1 Mpc and quasar 
lifetime t = 10 8 yr. 

using the ionizing equilibrium relation 



where S is the dark matter overdensity, F is the inten- 
sity of the radiation background, the factor of 1.7 is ap- 
proximately valid for the observed temperature-density re- 
lation (for isothermal gas the slope is 2 and the results 
are a lmost identical), and A is a constant of proportion- 
ality iHui fc Gnedhjri997t) . For the ionizing background F 
we take the local value at that position. We then map the 
neutral density field into redshift space by adding the pecu- 
liar velocity to the position at each grid cell and interpolate 
this field back to the grid. We ignore any additional ther- 
mal smoothing, since the grid cell size is already larger than 
the thermal broadening scale. Finally we compute the flux 
F — exp(— r). We vary the constant A until achieving the 
desired mean flux decrem ent, F, to match observations (see 
IMeiksin fc White) i2004h for a recent compilation). Figure 3] 
shows the resulting power spectra as a function of redshift. 
Comparison to observations reveals good qualitative agree- 
ment over the whole range of scales and redshifts. When 
comparing the spectra with and without background radi- 
ation fluctuations we always match first the mean flux for 
the two cases. 

For each line of sight we compute the flux distribution 
using the uniform background approximation and using the 
actual radiation field at each point. We then compute the 
power spectrum for each case, averaging over thousands lines 
of sight in each case. In figures|^|and|S|we present the results 
for these ratios for both (3 = 3.41 and ft = 2.58. In all cases 
we show the scenario where the quasar contribution to the 
background is 100%. If quasars contribute a fraction a, with 
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Figure 4. Flux power spectrum from the TPM simulation. The 
redshifts from top to bottom are z = 4.58, 3.88, 3.29, 2.82, 2.40 
and 2.05. The spectra agree qualitatively well with the observa- 
tions. 

the rest of the background uniform, then the effects in figures 
15161 are reduced roughly by this amount. 

There are several qualitative features in figures 15161 
worth pointing out. First is the fact that at high redshifts 
where the effect is significant the fluctuations in ionizing 
background enhance the power spectrum on large scales. A 
second important feature is the rapid decrease of the im- 
portance of fluctuations below z = 4. This is a consequence 
of the increase in the attenuation length r at t, which leads 
to more quasars contributing to a given position and thus 
reducing the fluctuations. 

For the S P SS po wer spectrum measurements 
^McDonald et all l2004al ) the relevant range is 
10" 3 s/km < k < 2 x 10" 2 s/km and 2.2 < z < 4.2. 
These observations do not go to sufficiently high redshift 
for these radiation background fluctuation effects to be 
clearly visible in the power spectra. There is no evidence 
for any enhancement of power on large scales at the highest 
redshifts (z — 4.2); note, how ever, that the statistica l 
errors are relatively large here (McDonal d et alJ 12004b). 
Some extreme models, like tcjso = 10 7 yr, (3 — 2.58 and 
r a tt = 53/i _1 Mpc produce a 20% effect at k = 10 -3 s/km 
and z = 3.88. These models are strongly disfavored by 
the data. Other models, like tq^o = 10 8 yr, /3 = 3.4 and 
r a tt — 106/i _1 Mpc, produce no effect at this redshift and 
are acceptable even if quasars dominate the background 
radiation field. The effects are even smalle r at lower redshift. 

Comparing our results to lCroftJ (120031) at z = 3, we note 
that the effect found there is somewhat larger and has the 
opposite sign. It is unclear what the source of the difference 
is, but one possibility is the hybrid hydrodynamic-dark mat- 
ter approach used in ICrofd (2003), which allows for only a 
qualitative estimate of the effect. Another difference in the 



Figure 5. Ratio of the power spectra from simulations with 
ionizing background fluctuations versus uniform background, for 
various quasar lifetimes and attenuation lengths. This is for bright 
end slope of the quasar luminosity function /3 = 3.41. We assume 
all of the ionizing background comes from quasars. 



V \ 1 - 2_ 

z=4.58 ^< 


- I I 1 1 

z=2.82 


10»yr/70Mpc/h 

10»yr/35Mpc/h 0.8 n 
10'yr/70Mpc/h 

. 10'yr/3SMpc/h 


10»yr/232Mpc/h 

10»yr/H6Mpc/h 

10'yr/232Mpc/h 

10 7 yr/116Mpc/h 


J g 1 1 1 — I — I I I 1 1 1[ 1 — I 

■\ ^001 0.01 

\ z=3.88 
■ sip 


1 III 1 I I Mill] 1 1 

0.001 0.01 

z=2.40 


10»yr/106Mpo/h 

10»yr/53Mpc/h 8n 

10'yr/106Mpc/h 

, 10'yr/53Mpc/h 


10»yr/340Mpc/h 

10»yr/170Mpc/h 

10 , yr/340Mpc/h 

. 10'yr/170Mpc/h 


/.I III 1 I I III 1 1 

0.001 0.01 

1.2- 

\^ z=3.29 


. I III 1 I I Mill] 1 1 

0.001 0.01 

z=2.05 


10"yr/162Mpc/h 

10"yr/81Mpc/h 0.8- 

10'yr/162Mpc/h 

10'yr/BlMpp/h 


10"yr/484Mpc/h 

p 10"yr/242Mpc/h 
10'yr/484Mpc/h 

10'yr/242Mpc/h 



0.001 0.01 0.001 0.01 

k[s/km] k[s/km] 



Figure 6. Same as figure fusing bright end slope of the quasar 
luminosity function fi = 2.58. 
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two treatments is the choice of quasar hosts . Whil e we use 
actual halos, overdensities were used in ICroftl(l2003ft . Finally 
there are also some differences in the adopted quasar lumi- 
n osity function . While redshift dependence was not explored 
in lCrofti (|5)03), our analysis suggests that a model that pro- 
duces a 20% effect at k = 10~ 3 s/km and z = 3 will lead to 
a large effect at z — 4 and will be excluded by the ab sence 
of any such effect in real data iMcDonald et al1l2004aT) . 

There could be additional sources of ionizing radiation 
such as Lyman break galaxies dSteidel et al.ll2003l) or rera- 
diation of photons from the absorbing sources in the forest 
l|Haardt fc Madaul Il9 9tf) . It is often assumed that due to 
their higher density these sources lead to uniform radiation 
field and so their effects on ionizing fluctuations need not 
be considered. However, these sources can still be clustered 
and if the large scale effects are coherent their contribution 
could be important. Here we address this by choosing all ha- 
los with masses between fO 11 /i -1 Mq and fO 12 hT 1 Mq as the 
hosts of LBGs. Their number density is 6 x f CT 3 /i 3 Mpc -3 , 
compa rable or somewhat larger than the density of known 
LBGs dSteidel et al.ll2003l) . We assign a luminosity propor- 
tional to halo mass and assume long lifetimes, comparable 
to the Hubble time at a given redshift. This model need not 
be correct, since LBGs could be starburst galaxies hosted by 
small halos r ather than sitting in th e most massive halos at 
that redshift (Somcrvillc ct al, 1200 ill , but we lack mass res- 
olution to resolve these small halos and test this alternative 
model. It is likely that the starburst model leads to more 
homogeneous radiation background. 

The rest of the model is the same as for quasars and in 
our calculations we assume all of ionizing background comes 
from these sources. In figure|2|we show the effect on the flux 
power spectrum as a function of redshift. We see that the 
redshift dependence is still present, but is weaker than for 
quasars. The effect is small except on largest scales, where 
it can reach 10% at redshifts of interest. Interestingly, over 
the relevant range of scales the effect leads to a suppres- 
sion of power, contrary to the quasar case. This could be a 
consequence of coherence between the galaxy field and Lya 
forest which comes into better light for these high density 
sources, unlike the case of quasars where shot noise was more 
important. 

In an actual analysis of the real data iMcDonald et all 
2004b) we tried adopting the templates of the worst case 
quasar scenario and allowing the fractional quasar contri- 
bution to vary smoothly with redshift (while the estimates 
of quasar contribution to the total background are uncer- 
tain, we do n ot expect the relative fr action to vary wildly 
with redshift jMeiksin fc White! 12003^. The d ata strongly 
disfavor these models iMcDonald et al.ll2004bT) . This could 
be because the worst case scenario quasar model adopted 
is too extreme, so other models in figures 15161 could still 
be acceptable as they have little or no effect on the flux 
power spectrum. Alternatively, the quasar effect should be 
combined with the contribution from galaxies, which leads 
to an opposite effect and so the two contributions partially 
cancel out. 
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Figure 7. Example of the effect on the flux power spectrum if the 
sources of ionizing radiation are galaxies with density comparable 
to that of LBGs. We find a suppression of power in this case. 



4 GALACTIC SUPERWINDS 

The effects of GSW on the Lya forest power spectrum have 
been explored before l|Croft et alJl2002al : iDesiacaues et alJ 
2003), where the wind simulation power spectra were com- 
pared to those with no winds. In this paper we perform a 
similar analysis using hydrodynamic simulations with and 
without GSW, but we also account for the allowed changes 
in the mean absorption, temperature-density relation, and 
filtering l ength. These are part o f our standard Lya forest 
analysis jMcDonald et al.ll2004bT) and our goal is to inves- 
tigate whether there are additional effects that go beyond 
the simple changes in the mean transmitted flux fraction, 
F, mean temperature, To, slope of the temperature-density 
relation, 7 — I, and filtering length, Uf- 

Our simulations (similar, but not the same as those 
in H2.21 are based on a fixed grid Eulerian cosmological 
hydrodynamic c ode with a TVD shock-capturing scheme 
JCen et alJl2004l) . They include the effects of cooling, heat- 
ing, star formation, and supernova (SN) feedback. The cos- 
mology adopted here is standard ACDM with Q m = 0.29, 
Q b = 0.047, H = 70km/s/Mpc, and cr 8 = 0.85. All of the 
simulations used here have a box size of lf/i _1 Mpc in co- 
moving units, 21 6 3 dark matt e r part icles, and gas followed 
on a 432 3 grid. See lCen et all (2004) for additional simula- 
tions and resolution tests. 

As described in lCen et alJ ll2004l) . there is no attempt 
to model the subgrid physics in these simulations, which 
would in principle determine how much of the SN energy 
is injected into the galaxies and how much of it can escape 
them. Instead it is assumed the wind energy (and mass) flux 
is proportional to the star formation rate M* and t he propor- 
tiona l ity constant is cali brated on observations jHeckmanl 
2003; Pcttini ct al. 2001). For energy output this implies 
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Figure 8. Ratio of the power spectra between strong wind sim- 
ulation and no wind simulation for z = 2,3, 4. Dashed line shows 
direct comparison, solid line the comparison after adjusting F, 
the temperature-density relation, and the filtering scale. 



Egsw = eGswc 2 M* with easw = 3 x 10 -6 . We assume 
a similar relation for the mass output. These are somewhat 
uncertain and we want to explore worst case scenarios, so 
we will analyze simulations where the outflow energy is in- 
creased by 5 to cgs w = 1-5 x 1 0~ 5 . 

As discussed in lCen et alJ (|2004T) . the wind simulations 
produce significant outflows which can heat up the IGM and 
pollute it with metals at a level consistent with observations. 
However, they appear to have little effect on the distribu- 
tion of neutral hydrogen. This is because GSW propagate 
preferentially in the directions of low density, so while they 
heat up the gas in voids, they do not affect the structure 
in higher density regions where most of the fluctuations in 
absor ption by neutral hydrogen are produced llTheuns et alJ 
120021) . 

To quantify these statements we explore explicitly the 
effect of GSW on the Lya forest flux power spectrum. This 
is shown in figure |H| for the more extreme value easw = 
1.5 x 10 -5 at 3 different redshifts. We do not show the results 
for the more realistic value cgsw — 3 x 10 -6 since the effects 
there are less significant. While we see that the effects for 
this more extreme case can be up to 10% over the range 
of observational interest, especially at z = 3, the effect is 
mostly due to the fact that the IGM temperature history 
and density dependence have changed. Without winds, the 
temperature parameters at z = (2, 3, 4) are T1.4 =(21699, 
17680, 13478) K and 7 - 1 =(0.454, 0.397, 0.497), while 
with maximum winds, they change to Ti.4 =(19563, 15012, 
14455) K and 7 — 1 =(0.479, 0.493, 0.477), where T 1A is the 
temperature at density 1.4 times the mean. 

To explore the effect of changes in the temperature- 
density relation, we recompute the power spectra after mod- 
ifying the temperatures in the simulations to impose differ- 



ent T\a and 7 — 1. We do this in the following way: First 
we verify that computing the neutral hydrogen density from 
the gas density using the ionizing equilibrium assumption 
leads to identical results as the simulation outputs them- 
selves. Then we compute T1.4 and 7 — 1 for a simu lation 
using a least absolute deviation fit JPress et al.lll992l) to all 
the points in the simulation with 1 < p/p < 2. This fit 
effectively uses the median rather than the mean tempera- 
ture, which is desirable because the mean temperature val- 
ues are heavily biased by rare regions with extremely high 
temperatures in simulations with GSW. For each gas cell we 
compute from its density the expected temperature for the 
original temperature-density relation in the simulation and 
for a modified temperature-density relation. We take the ra- 
tio between the two and multiply the actual value of the 
temperature in the cell by it, therefore preserving the scat- 
ter in the temperature-density relation present in the hydro- 
dynamic simulation. We then use the new temperatures to 
compute the neutral density using the ionizing equilibrium 
equation, and finally recompute the flux power spectrum. 

To see how much of the difference between simula- 
tions with and without winds is accounted for by sim- 
ple temperature-density relation changes, we recompute the 
power spectra in the simulations with GSW for a variety 
of modified temperature-density relations until we find one 
that matches best the power spectrum of the simulation 
without GSW. We find that the required changes in T1.4 
and 7 — 1 are in similar to the changes in the median-fitted 
values, but we did not restrict ourself to these in our analysis 
since it is not clear that they are the appropriate averages. 
For the purpose of the power spectrum analysis, finding the 
values that minimize the power spectrum differences is most 
relevant, since in our power spectrum analysis we always 
marginalize over T1.4 and 7 — 1. 

While we find good agreement between Pp in wind and 
no-wind simulations using this procedure where we only 
modify Ti.4 and 7 — 1, we can improve it further by not- 
ing that a change in the temperatures also changes the fil- 
tering length kp. This filtering length depends on the full 
temperature history of the IGM from reionization onwards. 
We do not try to make analytic predictions of what the 
filtering length should be given the temperature history of 
the IGM due to the difficulties in interpreting the scatter in 
temperat ure at a given density. I nstead we follow the pro- 
cedure in McDonald et al. (2004b), using hydro-PM simula- 
tions IGnedin fc Huill998T) to find the effect of changing the 
reionization redshift. We apply the change that minimizes 
the residual effect of the winds when combined with changes 
in Ti.4 and 7 — 1. Finally, we also allow small adjustments 
in the mean absorption level, F, which must be marginal- 
ized over in any cosmological fit. The final result is shown 
in figure |H] where we see that using the freedom of varying 
F, T1.4, 7 — 1, and kp cancels completely the effect of GSW 
in the simulations over the region of observational interest. 

The parameter variations we used to produce this 
agreement are not extreme. Beyond the change in median 
temperature-density relation estimated directly from the 
simulations, the modifications needed are as follows: We 
modified the filtering (kp) effect by 20% of the difference 
between two cases: the case when reionization heats the 
gas to 25000 K at z = 7 and the case when reioniza- 
tion heats the gas to 50000 K at z = 17 (we present the 
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change this way because we interpolated between these two 
models instead of actuall y running intermediate simulations 
jMcDonald et al J 12004th '). Note that we do not adjust k F 
independently at each redshift - the early-time thermal his- 
tory is modified and this effects each of the redshifts we 
show in a self-consistent way. At z — 2, we make the addi- 
tional small modifications AF = -0.001, AT1.4 = -800 K, 
and A7 — 1 = —0.03; at z = 3, we make the modifications 
AF = +0.006, AT1.4 = -1000 K, and A7 - 1 = -0.01; 
and finally, at z = 4, we make the changes AF — +0.0065, 
AT1.4 = +1000 K, and A7 - 1 = 0.05. All of these changes 
are small relative to the uncertainty in our measurements of 
these parameters (M cDonald et alj|200 4bj). 

While the simulation box is too small to be able to 
probe larger scales, the flatness of the final result on large 
scales suggests that these effects are only important on small 
scales and that these conclusions are unlikely to be modified 
if larger simulation boxes were used. Note that we find no 
effect both in the simulation with the expected GSW ampli- 
tude, as well as in the one where it was increased by a factor 
of 5. We only show the worse case scenario in figure |B] 

While these results are comforting, they should not be 
viewed as a conclusive proof that GSW have no effect on 
Lya forest flux power spectrum. For one thing, it is not 
clear that our simulations fare any better than other simu- 
lations in explaining the flux enha ncement observed in lines 
of sight very close to LBGs (lAdelberger et al]l2003l) . While 
no realistic simulations to date have been able to reproduce 
this result entirely, it is possible that it is a consequence of 
a statistical fluctuation, hosts of LBGs being less massive 
than usually assumed, galaxy proximity effect, selection ef- 
fects, or any combinati on of these. Moreover, the fact that 
iDesiacaues et all <l2003t) come close to explaining these ob- 
servations without affecting the flux power spectrum is com- 
forting. However, it is also possible that we are still missing 
some essential physics and that the winds are much more 
destructive of the IGM, although there are additional con- 
strai nts on this like the flux probability distribution statis- 
tics jMcDonald et all 120001) . Even in such a case it is not 
clear how much of the effect is degenerate with the parame- 
ters we already marginalized over, such as F, T1.4, 7 — 1, and 
kF- It is clear that more work is needed to exclude all pos- 
sible sources of contamination of the Lya forest flux power 
spectrum, even if the current modeling suggests that GSW 
are not problematic. 



5 CONCLUSIONS 

In this paper we have analyzed physical effects on the power 
spectrum of the Lya forest. We focus on effects that would 
not o therwise be part of the st andard simulation-based anal- 
ysis iMcDonald et al]l2004tj) . In section 2 we analyze the 
contribution from high density systems, especially those 
with damping wings, which are not simulated properly in 
the current generation of simulations. We show that their 
inclusion adds power on large scales and accounting for this 
leads to an increase in the deduce d slope of primordial fluc- 
tuations iMcDonald et a l. 2004b). The amplitude of the ef- 
fect can be estimated both from direct counting of these high 
density systems as well as using the templates in the power 
spectrum analysis itself, taking advantage of their specific 



form as a function of scale and redshift. At the moment the 
uncertainties in the for mer make the latter ap proach com- 
petitive or even better (iMcDonald et alJl2004bT) . suggesting 
that a better determination of the high density systems as a 
function of column density and redshift would be useful for 
accurate determination of this effect. But even within the 
current uncertainties it is clear that the effect is appreciable 
compared to the size of the statistical error given by the new 
data. 

In section 3 we analyze the effects of fluctuations in 
the ionizing background. We focus on quasars, which are 
believed to be the dominant source of ionizing background 
at these redshifts. These fluctuations lead to an enhance- 
ment of power on large scales which rapidly increases with 
redshift, because the photon attenuation length is signifi- 
cantly shorter at higher redshifts. This allows us to develop 
templates that can be used to identify the effect in the 
real data. The SPSS da ta show no evidence of this effect 
iMcDonald et alJl2004bl) . We also investigated a simple ex- 
ample that could correspond to galactic sources of radiation, 
matching the density of LBG galaxies. Contrary to naive ex- 
pectations we do not find that this more dense population of 
sources leads to a more homogeneous ionizing background, 
likely as a result of clustering of sources. The coherent na- 
ture of these sources and Lya forest leads to a suppression of 
power on large scales. The effect is relatively small on scales 
probed by SDSS and is somewhat redshift dependent, which 
allows it to be constrained from the data. 

We should caution that the current level of modeling 
remains rather simplistic and there are many possible sce- 
narios we did not explore, so we cannot conclusively rule 
out that fluctuations in ionizing background are unimpor- 
tant for the flux power spectrum statistics. To put things 
into perspective we note that the effect of damped systems, 
which is 20% at k ~ 10~ 3 s/km (figure |2J|, results in a 0.06 
change in the sl ope of the power spectru m, which is slightly 
more than 1-er iMcDonald et alJl2004bl) . In comparison to 
this the effects from ionizing fluctuations are below 10% at 
k ~ 10 _3 s/km. Thus these effects are unlikely to change the 
estimate of the slope by more than 1-er given the current 
uncertainties. In addition, they have stronger redshift de- 
pendence than damped systems, which makes it possible to 
identify and remove them. It is however important to under- 
stand them better if we want to extract all of the information 
present in the SDSS data, as the current constraints on the 
slope are limited by theoretical modeling uncertainties and 
could be improved further. 

In contrast to the effects of damped systems and ioniz- 
ing background fluctuations, which mostly change the power 
spectrum on large scales, galactic superwinds are predomi- 
nantly affecting small scales, where several other parameters 
are also important. In section 4 we have shown that galactic 
superwinds used in our simulations do not affect the Lya for- 
est flux power spectrum statistics after marginalization over 
the mean absorption level, temperature-density relation of 
the IGM, and its filtering scale. This is true even though 
the same GSW can heat up a significant fraction of the vol- 
ume and pollute t he IGM with me tals at a level consistent 
with observations jCen et alj l2004) . This again is not a con- 
clusive proof that winds are unimportant and there could 
still be physics missing in these simulations, but a defini- 
tive answer can only come by investigating it in simulations 
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which are sufficiently realistic and are able to satisfy all of 
the observational constraints. 

In summary, current models suggest that various phys- 
ical effects analyzed in this paper, which are not part of the 
current simulations, are not destructive for the power spec- 
trum statistics, even if they may weaken its predictive power. 
This conclusion is preliminary and more work is needed to 
investigate these and other effects before their impact is fully 
understood and corrected for in the analysis. One of the ef- 
fects we did not investigate here are fluctuations in the IGM 
temperature, which could be caused by pat chy helium reion- 
ization at these redshifts JZaldarriagal2002ft . This is unlikely 
to be a major effect, since temperature has only a minor ef- 
fect on the flux power spectrum statistics on large scales 
probed by SDSS, but this should be verified with explicit 
calculations. Progress on these issues will come from both 
theoretical and observational directions. Better theoretical 
modeling of the effects discussed in this paper, with simula- 
tions that include more realistic physics, will allow a better 
assessment of their impact. Improvements in observational 
tests of the Lya forest will be equally important. A few ex- 
amples of these that go beyond the flux power spectrum 
are of abundances of damped systems, correlations between 
Lya forest and quasars or galaxies and higher order cor- 
relations of Lya forest such as the bispectrum analysis or 
one-point distribution of the flux. All of these will allow to 
better constrain the contribution of these effects to the flux 
power spectrum and determine their effect on the cosmolog- 
ical constraints derived from Lya forest . 
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